Bi-level optimization-based infrared and visible light fusion method

ABSTRACT

The present invention proposes a bi-level optimization-based infrared and visible light fusion method, adopts a pair of infrared camera and visible light camera to acquire images, and relates to the construction of a bi-level paradigm infrared and visible light image fusion algorithm, which is an infrared and visible light fusion algorithm using mathematical modeling. Binocular cameras and NVIDIA TX2 are used to construct a high-performance computing platform and to construct a high-performance solving algorithm to obtain a high-quality infrared and visible light fusion image. The system is easy to construct, and the input data can be acquired by using stereo binocular infrared and visible light cameras respectively; the program is simple and easy to implement; and the fusion image is divided into an image domain and a gradient domain for fusion by means of mathematical modeling according to different imaging principles of infrared and visible light cameras.

TECHNICAL FIELD

The present invention belongs to the field of image processing and computer vision, adopts a pair of infrared camera and visible light camera to acquire images, and relates to the construction of a bi-level paradigm infrared and visible light image fusion algorithm, which is an infrared and visible light fusion algorithm using mathematical modeling.

BACKGROUND

The binocular stereo vision technology based on visible light band is developed to be relatively mature. Visible light imaging has rich contrast, color and shape information, so the matching information between binocular images can be obtained accurately and quickly so as to obtain scenario depth information. However, visible light band imaging has defects, and the imaging quality thereof is greatly reduced, for example, in strong light, fog rain, snow or night, which affects the matching precision. Therefore, the establishment of a color fusion system by using the complementarity of different band information sources is an effective way to produce more credible images in special environments. For example, a visible light band binocular camera and an infrared band binocular camera are used to constitute a multi-band stereo vision system, and the advantage of not being affected by fog, rain, snow and light of infrared imaging is used to make up for the deficiency of visible light band imaging so as to obtain more complete and precise fusion information.

The multi-modality image fusion technology is an image processing algorithm that uses the complementarity and redundancy between a plurality of images and adopts a specific algorithm or rule for fusion to obtain images with high credibility and better vision. Compared with the singularity of the mono-modality fusion image, multi-modality image fusion can better obtain the interactive information of images in different modalities, and gradually becomes an important means for disaster monitoring, unmanned driving, military monitoring and deep space exploration. The goal is to use the difference and complementarity of imaging of sensors with different modalities to extract the image information of each modality to the greatest extent and use source images of different modalities to fuse a composite image with abundant information and high fidelity. Therefore, the multi-modality image fusion will produce more comprehensive understanding and more accurate positioning of the image. In recent years, most of fusion methods are researched and designed based on the transform domain without considering the multi-scale detail information of images, resulting in the loss of details in the fused image, for example, the public patent CN208240087U, an infrared and visible light fusion system and image fusion device. Therefore, the present invention performs optimization solution after mathematical modeling of infrared and visible light, and realizes the enhancement of details and the removal of artifacts on the basis of retaining the effective information of infrared and visible light images.

SUMMARY

The present invention aims to overcome the defects of the prior art and provide a bi-level paradigm-based infrared and visible light image fusion algorithm. Through the design, bi-level paradigm-based mathematical modeling is carried out on infrared and visible light images, the images are solved in the image domain and the gradient domain, and finally, the real-time performance is achieved through GPU acceleration.

The present invention has the following specific technical solution: A bi-level optimization-based infrared and visible light fusion method, comprises the following steps:

1) Obtaining registered infrared and visible light images, and respectively calibrating each lens and jointly calibrating the respective systems of the visible light binocular camera and the infrared binocular camera;

1-1) Respectively calibrating the infrared camera and the visible light camera by the Zhangzhengyou calibration method to obtain internal parameters including focal length and principal point position and external parameters including rotation and translation of each camera;

1-2) Calculating the positional relationship of the same plane in the visible light image and the infrared image by using RT (rotation matrix and translation vector) obtained by joint calibration and the detected checker corners, and registering the visible light image to the infrared image by using a homography matrix;

2) Converting the color space of the visible light image from an RGB image to an HSV image, extracting the value information of the color image as the input of image fusion, and retaining the original hue and saturation;

2-1) In view of the problem that the visible light image has RGB three channels, converting the RGB color space to the HSV color space, wherein V is value, H is hue and S is saturation; and extracting the value information of the visible light image to be fused with the infrared image, and retaining the hue and saturation, wherein the specific conversion is shown as follows: R′=R/255 G′=G/255 B′=B/255 Cmax=max(R′,G′,B′) Cmin=min(R′,G′,B′) Δ=Cmax−Cmin V=Cmax

2-2) Extracting the V channel as the input of visible light, retaining H and S to the corresponding matrix, and retaining the color information for the subsequent color restoration after fusion.

3) Carrying out bi-level paradigm-based mathematical modeling on the input infrared image and the visible light image with the color space converted; the core of the idea is to establish two separate models, namely the upper-level sub-problem Leader and the lower-level sub-problem Follower to solve the same problem. First, the objective function and the problem constraint in optimization are assumed to be two participants in the game, wherein the objective function is regarded as the Leader and the problem constraint is the Follower. In this kind of competition, the next optimization of the Leader needs to consider the result of the Follower. The two competitors are defined as two composite minimization sub-problems:

${{{Leader}(l)}:{{\min\limits_{F}{{F - I}}_{2}^{2}} + {\gamma{{F - V}}_{2}^{2}}}}{{{Follower}(f)}:{{s.t.F} \in {{\underset{F}{argmin}{{{\nabla F} - {\nabla I}}}_{1}} + {\beta{{{\nabla F} - {\nabla V}}}_{1}}}}}$

wherein F represents the fused image, the infrared image and the visible light image are respectively represented by I, V, ∇ represents the operator for obtaining the gradient, and γ, β represent parameters of Leader and Follower respectively;

4) Solving the upper-level sub-problem to obtain the fusion result of the target in the image domain; the target result is obtained by solving the following formula:

$F_{l}^{k + 1} = {{\underset{F}{\arg\;\min}{{F - I}}_{2}^{2}} + {\gamma{{F - V}}_{2}^{2}}}$

wherein F₁ ^(k−1) represents the result of the upper-level problem; and the goal has a simple closed-form solution and thus is directly obtained by the closed-form solution shown as the following formula:

$F_{l}^{k + 1} = \frac{1 + {\gamma V}}{1 + \alpha}$

wherein F represents the fused image, I represents the infrared image, V represents the visible light image, γ represents the weight parameter of the visible light and infrared image, and α represents the weight variable parameter.

5) Solving the lower-level sub-problem to obtain the fusion result of the target in the gradient domain; it can be found that the problem is non-convex and non-smooth and difficult to solve directly, the alternating direction multiplier method is used to transform an unconstrained problem into a constrained problem by introducing auxiliary variables, and then the problem is solved under the framework;

The solving method is: first introducing two auxiliary variables u, w, wherein u=∇F−∇I, w=∇F−∇V; and through variable substitution, transforming to minimizing the following problem:

${\min\limits_{F,u,w}{u}_{1}} + {\beta{w}_{1}} + {\frac{\rho_{1}}{2}{{{\nabla F} - {\nabla I} - u - \frac{\lambda_{1}}{\rho_{1}}}}_{2}^{2}} + {\frac{\rho_{2}}{2}{{{\nabla F} - {\nabla V} - w - \frac{\lambda_{2}}{\rho_{2}}}}_{2}^{2}} + {\frac{\tau}{2}{{{\nabla F} - {\nabla{\overset{\hat{}}{F}}^{k}}}}_{2}^{2}}$

wherein V represents the gradient operator, λ₁, λ₂ are two multipliers, ρ₁, ρ₂ are parameters of a penalty term, and three sub-problems respectively about u, w, F are obtained through variable separation:

${\min\limits_{u}{u}_{1}} + {\frac{\rho_{1}}{2}{{u - {\nabla F} + {\nabla I} + \frac{\lambda_{1}}{\rho_{1}}}}_{2}^{2}}$ ${\min\limits_{w}{\beta{w}_{1}}} + {\frac{\rho_{2}}{2}{{w - {\nabla F} + {\nabla V} + \frac{\lambda_{2}}{\rho_{2}}}}_{2}^{2}}$ ${\min\limits_{F}{\frac{\rho_{1}}{2}{{{\nabla F} - {\nabla I} - u - \frac{\lambda_{1}}{\rho_{1}}}}_{2}^{2}}} + {\frac{\rho_{2}}{2}{{{\nabla F} - {\nabla V} - w - \frac{\lambda_{2}}{\rho_{2}}}}_{2}^{2}}$

5-1) For the update of F_(f) ^(k+1) of the lower-level problem, the closed-form solution of the formula is used to obtain:

$F_{f}^{k + 1} = \frac{{\rho_{1}{\nabla^{T}\left( {{\nabla I} + u + \frac{\lambda_{1}}{\rho_{1}}} \right)}} + {\rho_{2}{\nabla^{T}\left( {{\nabla V} + w + \frac{\lambda_{2}}{\rho_{2}}} \right)}} + {I*{\nabla{\overset{\hat{}}{F}}^{k}}}}{\rho_{1}{\nabla^{T}{\nabla{+ \rho_{2}}}}{\nabla^{T}{\nabla{+ I}}}*{\nabla^{T}\nabla}}$

5-2) The multipliers λ₁, λ₂ need updating after each iteration, and the specific update mode is as follows: λ₁ ^(k+1)=λ₁ ^(k)−ρ₁(∇F−∇V−w) λ₂ ^(k+1)=λ₂ ^(k)−ρ₂(∇F−∇I−u)

6) Obtaining two estimates F₁ ^(k+1), F_(f) ^(k+1) of the fusion result under different characteristics by solving the upper-level and lower-lever sub-problems, and to fuse the two components into an image F, linearly combining the two components, which is expressed as the following form: F ^(k+1)=α^(k) F _(l) ^(k+1)+(1−α^(k))F _(f) ^(k+1)

wherein F^(k+1) is the final result of each iteration, and α is a parameter weighing the two components and is manually selected based on experience.

7) Converting the color space: converting the fusion image back to the RGB image, and adding the hue and saturation previously retained;

Restoring to the RGB color space from the HSV color space by updating the V information saved into the fusion image in combination with the previously retained H and S, wherein the specific formulas are shown as follows:

C = V × SX = C × (1 − (H/60^(∘))mod 2 − 1) $m = {{V - {C\left( {R^{\prime},G^{\prime},B^{\prime}} \right)}} = \left\{ \begin{matrix} {\left( {C,X,0} \right)\ ,\ {{0{^\circ}} \leq H < {60{^\circ}}}} \\ {\left( {X,C,0} \right)\ ,\ {{60{^\circ}} \leq H < {120{^\circ}}}} \\ {\left( {0,C,X} \right)\ ,\ {{120{^\circ}} \leq H < {180{^\circ}}}} \\ {\left( {0,X,C} \right)\ ,\ {{180{^\circ}} \leq H < {240{^\circ}}}} \\ {\left( {X,0,C} \right)\ ,\ {{240{^\circ}} \leq H < {300{^\circ}}}} \\ {\left( {C,0,X} \right)\ ,\ {{300{^\circ}} \leq H < {360{^\circ}}}} \end{matrix} \right.}$ R′,G′,B′=((R′+m)×255,(G′+m)×255,(B′+m)×255)

wherein C is the product of the value and the saturation; and in is the difference of the value and C.

8) Enhancing the color: enhancing the color of the fusion image to generate a fusion image with higher resolution and contrast; and performing pixel-level image enhancement for the contrast of each pixel.

Performing color correction and enhancement on the restored image to generate a three-channel image that is consistent with observation and detection; and performing color enhancement on the R channel, G channel and B channel respectively to obtain the final fusion image. The specific formulas are shown as follows: R _(out)=(R _(in))^(1/gamma) R _(display)=(R _(in) ^((1/gamma)))^(gamma) G _(out)=(G _(in))^(1/gamma) G=(G _(in) ^((1/gamma)))^(gamma) B _(out)=(B _(in))^(1/gamma) B _(display)=(B _(in) ^((1/gamma)))^(gamma)

wherein gamma is the correction parameter, R_(in), G_(in) and B_(in) are the values of the three input channels R, G, and B respectively, R_(out), G_(out) and B_(out) are the intermediate parameters, and R_(display), G_(display) and B_(display) are the values of the three channels after enhancement.

The present invention has the following beneficial effects:

The present invention proposes a real-time fusion method using infrared and visible light binocular stereo cameras. Paradigm-based mathematical modeling is carried out on the input infrared and visible light images by means of mathematical modeling, and the images are solved and fused in the image domain and the gradient domain respectively, effectively reducing the interference of artifacts and fusing a highly reliable image. The present invention has the following characteristics:

(1) The system is easy to construct, and the input data can be acquired by using stereo binocular cameras;

(2) The program is simple and easy to implement;

(3) The fusion is divided into two parts and specifically solved by means of mathematical modeling;

(4) The structure is complete, multi-thread operation can be performed, and the program is robust;

(5) The detail images are used to perform significant enhancement and differentiation to improve the generalization ability of the algorithm.

DESCRIPTION OF DRAWINGS

FIG. 1 is a flow chart of a visible light and infrared fusion algorithm.

FIG. 2 is a final fusion image.

DETAILED DESCRIPTION

The present invention proposes a method for real-time image fusion by an infrared camera and α visible light camera, and will be described in detail below in combination with drawings and embodiments.

The binocular stereo cameras are placed on a fixed platform, the image resolution of the experiment cameras is 1280×720, and the field of view is 45.4°. To ensure real-time performance, NVIDIA TX2 is used for calculation. On this basis, a real-time infrared and visible light fusion method is designed, and the method comprises the following steps:

1) Obtaining registered infrared and visible light images;

1-1) Respectively calibrating each lens and jointly calibrating the respective systems of the visible light binocular camera and the infrared binocular camera;

1-2) Respectively calibrating the infrared camera and the visible light camera by the Zhangzhengyou calibration method to obtain internal parameters such as focal length and principal point position and external parameters such as rotation and translation of each camera.

1-3) Calculating the positional relationship of the same plane in the visible light image and the infrared image by using RT obtained by joint calibration and the detected checker corners, and registering the visible light image to the infrared image by using a homography matrix.

2) Converting the color space of the image

2-1) In view of the problem that the visible light image has RGB three channels, converting the RGB color space to the HSV color space, extracting the V (value) information of the visible light image to be fused with the infrared image, and retaining H (hue) and S (saturation), wherein the specific conversion is shown as follows: R′=R/255 G′=G/255 B′=B/255 Cmax=max(R′,G′,B′) Cmin=min(R′,G′,B′) Δ=Cmax−Cmin V=Cmax

2-2) Extracting the V (value) channel as the input of visible light, retaining H (hue) and S (saturation) to the corresponding matrix, and retaining the color information for the subsequent color restoration after fusion.

3) Carrying out bi-level paradigm-based mathematical modeling for the input infrared image and the visible light image with the color space converted. The core of the idea is to establish two separate models, namely the upper-level sub-problem Leader and the lower-level sub-problem Follower to solve the same problem. First, the objective function and the problem constraint in optimization are assumed to be two participants in the game, wherein the objective function is regarded as the Leader and the problem constraint is the Follower. In this kind of competition, the next optimization of the Leader needs to consider the result of the Follower. The two competitors are defined as two composite minimization sub-problems:

${{{Leader}(l)}:{{\min\limits_{F}{{F - I}}_{2}^{2}} + {\gamma{{F - V}}_{2}^{2}}}}{{{Follower}(f)}:{{s.t.F} \in {{\underset{F}{argmin}{{{\nabla F} - {\nabla I}}}_{1}} + {\beta{{{\nabla F} - {\nabla V}}}_{1}}}}}$

wherein F represents the fused image, and the infrared image and the visible light image are respectively represented by I, V, ∇ represents the operator for obtaining the gradient, and γ, β represent parameters of Leader and Follower respectively.

4) Solving the upper-level sub-problem; the target result can be obtained by solving the following formula:

$F_{1}^{k + 1} = {{\underset{F}{argmin}{{F - I}}_{2}^{2}} + {\gamma{{F - V}}_{2}^{2}}}$

wherein F₁ ^(k−1) represents the result of the upper-level problem. The goal has a simple closed-form solution and thus can be directly obtained by the closed-form solution shown as the following formula:

$F_{1}^{k + 1} = \frac{1 + {\gamma V}}{1 + \alpha}$

5) Solving the lower-level sub-problem; it can be found that the problem is non-convex and non-smooth and difficult to solve directly and is solved by the alternating direction multiplier method widely used. Therefore, it is necessary to transform an unconstrained problem into a constrained problem by introducing auxiliary variables, and then the problem is solved under the framework. Specifically, first introducing two auxiliary variables u, w, wherein u=∇F−∇I, w=∇F−∇V. Through variable substitution, transforming to minimizing the following problem:

${\min\limits_{F,u,w}{u}_{1}} + {\beta{w}_{1}} + {\frac{\rho_{1}}{2}{{{\nabla F} - {\nabla I} - u - \frac{\lambda_{1}}{\rho_{1}}}}_{2}^{2}} + {\frac{\rho_{2}}{2}{{{\nabla F} - {\nabla V} - w - \frac{\lambda_{2}}{\rho_{2}}}}_{2}^{2}} + {\frac{\tau}{2}{{{\nabla F} - {\nabla{\overset{\hat{}}{F}}^{k}}}}_{2}^{2}}$

wherein ∇ represents the gradient operator, λ₁, λ₂ are two multipliers, ρ₁, ρ₂ are parameters of a penalty term, and three sub-problems respectively about u, w, F are obtained through variable separation:

${{\min\limits_{u}{u}_{1}} + {\frac{\rho_{1}}{2}{{u - {\nabla F} + {\nabla I} + \frac{\lambda_{1}}{\rho_{1}}}}_{2}^{2}}}{{\min\limits_{w}{\beta{w}_{1}}} + {\frac{\rho_{2}}{2}{{w - {\nabla F} + {\nabla V} + \frac{\lambda_{2}}{\rho_{2}}}}_{2}^{2}}}{{\min\limits_{F}{\frac{\rho_{1}}{2}{{{\nabla F} - {\nabla I} - u - \frac{\lambda_{1}}{\rho_{1}}}}_{2}^{2}}} + {\frac{\rho_{2}}{2}{{{\nabla F} - {\nabla V} - w - \frac{\lambda_{2}}{\rho_{2}}}}_{2}^{2}}}$

5-1) For the update of F_(f) ^(k+1) of the lower-level problem, the closed-form solution of the formula is used to obtain:

$F_{f}^{k + 1} = \frac{{\rho_{1}{\nabla^{T}\left( {{\nabla I} + u + \frac{\lambda_{1}}{\rho_{1}}} \right)}} + {\rho_{2}{\nabla^{T}\left( {{\nabla V} + w + \frac{\lambda_{2}}{\rho_{2}}} \right)}} + {I*{\nabla{\overset{\hat{}}{F}}^{k}}}}{\rho_{1}{\nabla^{T}{\nabla{+ \rho_{2}}}}{\nabla^{T}{\nabla{+ I}}}*{\nabla^{T}\nabla}}$

5-2) The multipliers λ₁, λ₂ need updating after each iteration, and the specific update mode is as follows: λ₁ ^(k+1)=λ₁ ^(k)−ρ₁(∇F−∇V−w) λ₂ ^(k+1)=λ₂ ^(k)−ρ₂(∇F−∇I−u)

6) Obtaining two estimates F₁ ^(k+1), F_(f) ^(k+1) of the fusion result under different characteristics by solving the upper-level and lower-lever sub-problems, and to fuse the two components into an image F, linearly combining the two components, which is expressed as the following form: F ^(k+1)=α^(k) F _(l) ^(k+1)+(1−α^(k))F _(f) ^(k+1)

wherein F is the final result of each iteration, and α is a parameter weighing the two components. The parameter needs to be manually selected based on experience, and is selected as 0.5 herein.

7-1) Restoring to the RGB color space from the HSV color space by updating the V (value) information saved into the fusion image in combination with the previously retained H (hue) and S (saturation), wherein the specific formulas are shown as follows:

${C = {V \times S}}{X = {C \times \left( {1 - {❘{{\left( {{H/60}{^\circ}} \right){mod}2} - 1}❘}} \right)}}{m = {V - C}}{\left( {R^{\prime},G^{\prime},B^{\prime}} \right) = \left\{ \begin{matrix} \left( {C,X,0} \right) & {{0{^\circ}} \leq H < {60{^\circ}}} \\ \left( {X,C,0} \right) & {{60{^\circ}} \leq H < {120{^\circ}}} \\ \left( {0,C,X} \right) & {{120{^\circ}} \leq H < {180{^\circ}}} \\ \left( {0,X,C} \right) & {{180{^\circ}} \leq H < {240{^\circ}}} \\ \left( {X,0,C} \right) & {{240{^\circ}} \leq H < {300{^\circ}}} \\ \left( {C,0,X} \right) & {{300{^\circ}} \leq H < {360{^\circ}}} \end{matrix} \right.}$ R′,G′,B′=((R′+m)×255,(G′+m)×255,(B′+m)×255)

wherein C is the product of the value and the saturation; and in is the difference of the value and C.

7-2) performing color correction and enhancement on the image restored in step 7-1 to generate a three-channel image that is consistent with observation and detection; and performing color enhancement on the R channel, G channel and B channel respectively, wherein the specific formulas are shown as follows: R _(out)=(R _(in))^(1/gamma) R _(display)=(R _(in) ^((1/gamma)))^(gamma) G _(out)=(G _(in))^(1/gamma) G=(G _(in) ^((1/gamma)))^(gamma) B _(out)=(B _(in))^(1/gamma) B _(display)=(B _(in) ^((1/gamma)))^(gamma)

wherein gamma is the correction parameter, R_(in), G_(in) and B_(in) are the values of the three input channels R, G, and B respectively, R_(out), G_(out) and B_(out) are the intermediate parameters, and R_(display), G_(display) and B_(display) are the values of the three channels after enhancement. 

The invention claimed is:
 1. A bi-level optimization-based infrared and visible light fusion method, wherein the method comprises the following steps: 1) obtaining registered infrared and visible light images, and respectively calibrating each lens and jointly calibrating the respective systems of the visible light binocular camera and the infrared binocular camera; 1-1) respectively calibrating the infrared camera and the visible light camera by the Zhangzhengyou calibration method to obtain internal parameters including focal length and principal point position and external parameters including rotation and translation of each camera; 1-2) calculating the positional relationship of the same plane in the visible light image and the infrared image by using RT obtained by joint calibration and the detected checker corners, and registering the visible light image to the infrared image by using a homography matrix; 2) converting the color space of the visible light image from an RGB image to an HSV image, extracting the value information of the color image as the input of image fusion, and retaining the original hue and saturation; 3) carrying out bi-level paradigm-based mathematical modeling on the input infrared image and the visible light image with the color space converted; establishing two separate models, namely the upper-level sub-problem Leader and the lower-level sub-problem Follower to solve the same problem. ${{{Leader}(l)}:{{\min\limits_{F}{{F - I}}_{2}^{2}} + {\gamma{{F - V}}_{2}^{2}}}}{{{Follower}(f)}:{{s.t.F} \in {{\underset{F}{argmin}{{{\nabla F} - {\nabla I}}}_{1}} + {\beta{{{\nabla F} - {\nabla V}}}_{1}}}}}$ wherein F represents the fused image, the infrared image and the visible light image are respectively represented by I, V, ∇ represents the operator for obtaining the gradient, and γ, β represent parameters of Leader and Follower respectively; 4) solving the upper-level sub-problem to obtain the fusion result of the target in the image domain; the target result is obtained by solving the following formula: $F_{1}^{k + 1} = {{\underset{F}{argmin}{{F - I}}_{2}^{2}} + {\gamma{{F - V}}_{2}^{2}}}$ wherein F₁ ^(k+1) represents the result of the upper-level problem; and the goal has a simple closed-form solution and thus is directly obtained by the closed-form solution shown as the following formula: $F_{1}^{k + 1} = \frac{1 + {\gamma V}}{1 + \alpha}$ wherein F represents the fused image, I represents the infrared image, V represents the visible light image, γ represents the weight parameter of the visible light and infrared image, and α represents the weight variable parameter; 5) solving the lower-level sub-problem to obtain the fusion result of the target in the gradient domain; and using the alternating direction multiplier method to transform an unconstrained problem into a constrained problem by introducing auxiliary variables, and then solving the problem under the framework; 6) obtaining two estimates F₁ ^(k+1), F_(f) ^(k+1) of the fusion result under different characteristics by solving the upper-level and lower-lever sub-problems, and to fuse the two components into an image F, linearly combining the two components, which is expressed as the following form: F ^(k+1)=α^(k) F _(l) ^(k+1)+(1−α^(k))F _(f) ^(k+1) wherein F^(k+1) is the final result of each iteration and is manually selected based on experience; 7) converting the color space: converting the fusion image back to the RGB image, and adding the hue and saturation previously retained; restoring to the RGB color space from the HSV color space by updating the V information saved into the fusion image in combination with the previously retained H and S; 8) enhancing the color: enhancing the color of the fusion image to generate a fusion image with higher resolution and contrast; and performing pixel-level image enhancement for the contrast of each pixel; performing color correction and enhancement on the restored image to generate a three-channel image that is consistent with observation and detection; and performing color enhancement on the R channel, G channel and B channel respectively to obtain the final fusion image.
 2. The bi-level optimization-based infrared and visible light fusion method according to claim 1, wherein the color space conversion of the visible light image in step 2) comprises: 2-1) converting the RGB color space to the HSV color space, wherein V is value, H is hue and S is saturation; and extracting the value information of the visible light image to be fused with the infrared image, and retaining the hue and saturation, wherein the specific conversion is shown as follows: R′=R/255 G′=G/255 B′=B/255 Cmax=max(R′,G′,B′) Cmin=min(R′,G′,B′) Δ=Cmax−Cmin V=Cmax 2-2) extracting the V channel as the input of visible light, retaining H and S to the corresponding matrix, and retaining the color information for the subsequent color restoration after fusion.
 3. The bi-level optimization-based infrared and visible light fusion method according to claim 1, wherein the solving method in step 5): first introducing two auxiliary variables u, w, wherein u=∇F−∇I, w=∇F−∇V; and through variable substitution, transforming to minimizing the following problem: ${\min\limits_{F,u,w}{u}_{1}} + {\beta{w}_{1}} + {\frac{\rho_{1}}{2}{{{\nabla F} - {\nabla I} - u - \frac{\lambda_{1}}{\rho_{1}}}}_{2}^{2}} + {\frac{\rho_{2}}{2}{{{\nabla F} - {\nabla V} - w - \frac{\lambda_{2}}{\rho_{2}}}}_{2}^{2}} + {\frac{\tau}{2}{{{\nabla F} - {\nabla{\overset{\hat{}}{F}}^{k}}}}_{2}^{2}}$ wherein ∇ represents the gradient operator, λ₁, λ₂ are two multipliers, ρ₁, ρ₂ are parameters of a penalty term, and three sub-problems respectively about u, w, F are obtained through variable separation: ${{\min\limits_{u}{u}_{1}} + {\frac{\rho_{1}}{2}{{u - {\nabla F} + {\nabla I} + \frac{\lambda_{1}}{\rho_{1}}}}_{2}^{2}}}{{\min\limits_{w}{\beta{w}_{1}}} + {\frac{\rho_{2}}{2}{{w - {\nabla F} + {\nabla V} + \frac{\lambda_{2}}{\rho_{2}}}}_{2}^{2}}}{{\min\limits_{F}{\frac{\rho_{1}}{2}{{{\nabla F} - {\nabla I} - u - \frac{\lambda_{1}}{\rho_{1}}}}_{2}^{2}}} + {\frac{\rho_{2}}{2}{{{\nabla F} - {\nabla V} - w - \frac{\lambda_{2}}{\rho_{2}}}}_{2}^{2}}}$ 5-1) for the update of F_(f) ^(k+1) of the lower-level problem, the closed-form solution of the formula is used to obtain: $F_{f}^{k + 1} = \frac{{\rho_{1}{\nabla^{T}\left( {{\nabla I} + u + \frac{\lambda_{1}}{\rho_{1}}} \right)}} + {\rho_{2}{\nabla^{T}\left( {{\nabla V} + w + \frac{\lambda_{2}}{\rho_{2}}} \right)}} + {I*{\nabla{\overset{\hat{}}{F}}^{k}}}}{\rho_{1}{\nabla^{T}{\nabla{+ \rho_{2}}}}{\nabla^{T}{\nabla{+ I}}}*{\nabla^{T}\nabla}}$ 5-2) the multipliers λ₁, λ₂ need updating after each iteration, and the specific update mode is as follows: λ₁ ^(k+1)=λ₁ ^(k)−ρ₁(∇F−∇V−w) λ₂ ^(k+1)=λ₂ ^(k)−ρ₂(∇F−∇I−u).
 4. The bi-level optimization-based infrared and visible light fusion method according to claim 1, wherein the specific formulas for color space conversion in step 7) are shown as follows: ${C = {V \times S}}{X = {C \times \left( {1 - {❘{{\left( {{H/60}{^\circ}} \right){mod}2} - 1}❘}} \right)}}{m = {V - C}}{\left( {R^{\prime},G^{\prime},B^{\prime}} \right) = \left\{ \begin{matrix} \left( {C,X,0} \right) & {{0{^\circ}} \leq H < {60{^\circ}}} \\ \left( {X,C,0} \right) & {{60{^\circ}} \leq H < {120{^\circ}}} \\ \left( {0,C,X} \right) & {{120{^\circ}} \leq H < {180{^\circ}}} \\ \left( {0,X,C} \right) & {{180{^\circ}} \leq H < {240{^\circ}}} \\ \left( {X,0,C} \right) & {{240{^\circ}} \leq H < {300{^\circ}}} \\ \left( {C,0,X} \right) & {{300{^\circ}} \leq H < {360{^\circ}}} \end{matrix} \right.}$ R′,G′,B′=((R′+m)×255,(G′+m)×255,(B′+m)×255) wherein C is the product of the value and the saturation; and in is the difference of the value and C.
 5. The bi-level optimization-based infrared and visible light fusion method according to claim 1, wherein the specific formulas for color enhancement in step 8) are shown as follows: R _(out)=(R _(in))^(1/gamma) R _(display)=(R _(in) ^((1/gamma)))^(gamma) G _(out)=(G _(in))^(1/gamma) G=(G _(in) ^((1/gamma)))gamma B _(out)=(B _(in))¹/gamma B _(display)=(B _(in) ^((1/gamma)))^(gamma) wherein gamma is the correction parameter, R_(in), G_(in) and B_(in) are the values of the three input channels R, G, and B respectively, R_(out), G_(out) and B_(out) are the intermediate parameters, and R_(display), G_(display) and B_(display) are the values of the three channels after enhancement. 